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ASYMPTOTIC BEHAVIOUR OF VISCOELASTIC COMPOSITES 
WITH ALMOST PERIODIC MICROSTRUCTURES 

JEAN LOUIS WOUKENG 


Abstract. In this paper we study the acoustic properties of porous media 
saturated by an incompressible viscoelastic fluid. The model considered here 
consists of a linear deformable porous skeleton having memory that is sur¬ 
rounded by a viscoelastic Oldroyd fluid. Assuming the microstructures to be 
almost periodically distributed and under the almost periodicity hypothesis 
on the coefficients of the governing equations, we determine the macroscopic 
equivalent medium. To achieve our goal, we use some very recent tools about 
the sigma convergence of convolution sequences. 


1. Introduction and the model 

A number of important engineering materials such as polymers store and dis¬ 
sipate simultaneously mechanical energy when subjected to external forces. This 
form of response as a combination of both liquid-like and solid-like features is termed 
viscoelasticity. The constitutive equations for the viscoelastic behaviour include the 
elastic deformation and viscous flow as special cases, and at the same time provide 
responses which are characteristic of the combined behaviour. In order to deal 
with viscoelastic fluid motion, researchers have focused on the use of constitutive 
equations that capture important qualitative features describing the relationships 
between the kinematic, mechanical and thermal fields equations and hence permit 
the formulation of well-posed problems that are simplistic relative the formulations 
that are supposed to characterize real materials. In that direction, many models 
have been proposed and developed. We refer e.g. to [D 111 [71 uni m m m [23 ISO] , 
just to cite a few. It is worth noting that in all the previous works, the models 
involved did not deal with skeletons with memory. 

In this work, we mean to study the acoustic properties of porous saturated media. 
In this respect we consider the problem of viscoelastic Oldroyd fluid motion through 
linear deformable porous media with memory. The novelty is threefold. First, 
the deformable skeleton has memory. Second, instead of considering usual Stokes 
equations flow, we rather consider the Oldroyd class of incompressible viscoelastic 
fluids. Third, the pores are distributed in an almost periodic fashion, generalizing 
the usual equidistribution (also known as periodic distribution) always considered 
in the earlier works. To be precise, as stated before, the problem will be discussed in 
terms of two particular types of ideal materials: the elastic skeleton (with memory) 
and the viscoelastic fluid, and the constitutive equations will therefore be described 
on the microscale as follows. 
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1.1. Description of the domain. Let Y = (0,1)^ (iV = 2 or 3 being fixed once 
for all here and in the sequel) be the reference cell, and let Yi and V 2 be two open 
disjoint subsets of V representing the local structure of the skeleton (which is a 
linear elastic material) and the local structure of the fluid, respectively. We assume 
that Yi C Y, Y = Yi U Y 2 , Y 2 is connected and that the boundary dYi of Yi is 
Lipschitz continuous. Moreover we assume that Yi and Y 2 have positive Lebesgue 
measure. Let S C be an infinite subset of Z^ and let 

Gi = Uk^sik + Yi), G2 =R^\Gi. 

Then Gj {j = 1, 2) is an open subset of K.^ and moreover G 2 is connected. Since 
the cells {k + Y)kes are pairwise disjoint, the characteristic function Xi of the set 
Gi in satisfies Xi = '^kes Xk+Yi or more precisely, 

Xl= YI 

fcez'v 

where 9 is the characteristic function of S in Z^. The function 6 is the distribution 
function of the pores (that is, it characterizes the porous cells; remind that not all 
the cells are porous). The special case 9{k) = 1 for all k € (i.e. S = 

leads to the equidistribution of fissured cells, commonly known in the literature as 
the periodic porous medium. Here we assume the function 9 to be almost periodic 
(which corresponds to the almost periodic distribution of the pores, that is, of the 
porous cells), i.e., the translates 9{- Y h) {h € Z^) form a relatively compact set in 
= L°°(Z^,A), A being the counting measure. Here the norm in 
is given by ||u||^oo( 2 'V) = sup^^ggN \u{k)\ for u € As seen in the remark 

below, the almost periodicity assumption encompasses all kinds of periodicity. 

Remark 1.1. 1) We recall that the space £°^p{'L^) in which the function 9 lies, is 
the smallest closed subspace of £°°{'E^) that contains all sequences of the form 

C : k 1 -^ exp(t^ ■ k), p, G 

Therefore it contains in particular all the periodic sequences defined on the lattice 
Z^ . 2) The function 9 being the characteristic function of a subset of Z^, assuming 
9 to be almost periodic in the above sense yields 9 is actually periodic. Indeed we 
find that there is a p G Z^ such that 

sup \9{k + p) - 9{k)\ < i. 
feez" ^ 

Since 9 assumes only 0 or 1, it follows that \9{k + p) — 9{k)\ takes on either 0 or 1. 
We deduce at once that 9{k + p) — 9{k) = 0 for all k G Z^ , so that 9 is periodic. 
Here, what we gain is that 9 may assume any kind of periodicity. 


Bearing this in mind, let be an open bounded Lipschitz domain in and 
let e > 0 be a small parameter. We define the almost periodic porous saturated 
medium as follows: (j = 1,2). The volumes and HI are respectively 

filled by the solid and the fluid, and H| is connected. We have H = Hf U rf 2 U H| 
(disjoint union) where rf 2 = n cIHI fl H is the part of the interface of Hf with 
HI lying in H. Let Vj denote the unit outward normal on 9H^; then vi = —V 2 on 
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1.2. The constitutive equations. In the solid. The elastic constitutive law cou¬ 
pled with the momentum balance give rise to 

d'^u r /■* 

AgVUe -I- / Al{x,t — T)\7Ue{x,T)dT 

Jo 


Pi- 


dt^ 


— div 


= plf in ni X (0,r) (1.1) 


where is the solid displacement, Aj the elastic tensors of the material constituting 
the skeleton and p\ the density of that material. 

In the fluid. The viscoelastic constitutive law associated to the momentum bal¬ 
ance and the incompressibility give rise to 


j- 


BqVVs+ f Bl{x,t - T)WVs{x,T)dT 

Jo 


+ Vpe = P29 in ^^2 X (0: T) 


( 1 . 2 ) 

divtig = 0 in fl 2 X (0,T) (1-3) 

where Vg denotes the velocity of the fluid, the dynamical viscosity tensors of 
the fluid, pI th® density of the fluid and the pressure of the fluid. 

At the fluid-solid interface. The adherence condition 


V , = -^011 rf 2 x (0, T ) 


and the normal stress continuity 

rt 


-I- J Al{x,t - T)VUsdT^ • Ui 

-I- J i?i(x,t — r)Vt)edr^ • ui on rf 2 X (0,T) 


(1.4) 

(1.5) 


= -Pel + 


hold. 


Finally we supplement Eq. (1III])-(II31) by the boundary conditions 

«£ = 0 on (cinf n (?n) x (O, r), = O on bl dfl) X (0, T) 


( 1 . 6 ) 


and the initial homogeneous conditions 

du 

Ug{x, 0) = 0) = 0 in flf and Ue(x, 0) = 0 in (1.7) 

Here, Aq{x) = Ao{x/e) and Hf(x,t) = Ai{xfe,tfe), Bq{x) = Bo{x/e) and 
Bl{x,t) = Bxixjeflje) ((x,t) e Q = H x (0,T), T > 0 a given real number), 
P]{x) = Pjixle) (j = 1,2) with: 

(Al) Aq, Bq e L°°(IR^)^^ and Ai,Bi G L°°(R^+^)^^ being symmetric matrices 
satisfying the following assumption: there exists a > 0 such that 

^ 0 ^ • ^ > a and Bof ■ f > a a.e. in 11 and for all ^ G R^. 

(A2) The density function pj (j = 1,2) lies in {Ry ) and satisfy A“^ < p^ (y) < 
A a.e. y G R^, for some positive A. 

(A3) The functions / and g belong to 

In view of the assumptions on Ai,Bi and pj (1 = 0,1, j = 1, 2), the functions Aq 
( resp. Bq), Af (resp. Bf) and pj are well defined and belong to (resp. 

L“(fl)^'), L°°(Q)^' (resp. L°°(Q)^') and L“(fl). 

We stress that we may view our model as a physically correct one in the sense 
that if in Equations (HD and (11.31) (and in the interface conditions (11.41) and (11.51) 1 
we replace the gradients Vitg and Vv^ by the strains e{uflj = i (Vug -I- S/^u^) 
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and e{ve) = \ (Vm^ + respectively, then thanks to Korn’s inequality [HI p. 

31], the mathematical analysis does not change, and the model becomes therefore 
physically correct. Thus the use of the gradient instead of the strain is just to 
simplify the exposition of the results. Also, the fact that the coefficients Aq and 
Bq do not depend on the fast time variable is just to avoid tricky and useless 
complications. Indeed assuming for example Aq dependent on t/e, we will be led 
to an estimation containing a term of the order 

The model in (fTTt-ifra is a generalization of those in 11017]. It seems to 
be more accurate to study the acoustic properties of microstructured materials 
such as glass wool, foam, etc. Let us precise that here, our considerations differ 
from those performed in all the previous work. Indeed, first, in the framework of 
deterministic homogenization theory, all the problems considered so far have been 
studied in periodic porous media. Here we consider an environment that is almost 
periodically perforated, generalizing the results of the earlier works. Second, the 
problems solved so far did not involve any memory effect at the microscale level; 
see e.g., [Ii|2Ql|35]. Here we consider a problem modeling a phenomenon with a 
memory effect at the fast time-variable. From the physical point of view, it means 
that the memory effects arising by meeting an obstacle decay in the surrounding of 
the next obstacle. So considered, the problem is new and the obtained results are 
of physical interest. Finally, though the problem considered here is linear, we do 
not use the Laplace transform to achieve our goal. We rather use a direct method, 
the two-scale convolution (in the general framework of ^-convergence)^ leaning on 
some recent results about the convergence of sequences defined by the convolution; 
see e.g.. Theorem 14.51 Let us precise that the periodic counterpart of Theorem 
14.51 is already known and has been proved by Visintin [HI Proposition 2.13] by 
using the periodic unfolding method. A more general version of that result in the 
framework of S-convergence has been for the first time stated and proved in m, 
but with the restricted assumption that the open set fl C is bounded. Still in 
the S-convergence framework, a more general proof (in any open set H C M.^) has 
been given in the almost periodic setting in a very recent work [32] . The use of 
this method instead of the Laplace transform method is motivated by at least two 
facts: 1) It enables one to tackle both linear and nonlinear homogenization problems 
involving convolution both in fast space and time variables (see EZllH]), which is 
not the case for the latter method; 2) As observed in numerical experiments by 
the authors of [35] (see the Section 5 therein), memory effects make numerical 
inverse Laplace transform very complicated''. 

Our aim is to pass to the limit as the size e of the microstructures shrinks to 0 
in the above system, and of course when the elasticity and viscosity tensors, the 
density of both fluid and skeleton, are subjected to be almost periodic, together 
with the assumption that the medium is almost periodically perforated. Based on 
the above considerations, we may therefore say that our results are new. 

The work is framed as follows. Section 2 deals with the statement and the proof 
of the existence result for problem dlllll-IIlTl). We also derive therein some a priori 
estimates and a compactness result. In section 3, we give some fundamental facts 
about almost periodicity. The sigma convergence concept with some related convo¬ 
lution results are discussed in Section 4. In Section 5, we derive the homogenization 
result leading to the properties of the equivalent medium as stated earlier. 
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Throughout the work, vector spaces are assumed to be real vector spaces, and 
scalar functions are assumed to take real values. We shall always assume that the 
numerical space (integer > 1) and its open sets are each equipped with the 
Lebesgue measure dx = dxi...dxN- 


2 . Existence and compactness results 

We begin by giving the functional setting of our problem. We define the function 
spaces 

: V e with the norm ||v||^, = ||Vu||^ 2 (opiv 2 , j = 1,2; 

W = {v ^ : divti = 0 in with usual IJq (n)^-norm; 

H = the closure of W for the norm induced by the inner product ((•,•)), 
where ((•,•)) denotes the weighted inner product 

((«,-*;))=/’ plu-vdx+ [ p^ 2 U-vdxioTU,v€L'^in)^. 

Jqi 

The induced norm on H is equivalent to the L^-norm and H preserves the diver¬ 
gence free property in ill- We denote by ((-,•)) the duality pairing between W 
and W generated from the above inner product ((-, •)). This being so, we have the 
evolution triple W ^ H ^ W' where we identify with its topological dual Id'. 

In order to solve (lOl-dTTTl). we need an auxiliary equivalent formulation of the 
said-problem. Prior to that, let us introduce the auxiliary functions 

w,(t) = [ v^{s)ds, = x\ue + xlwe, p" = xIpI + X 2 P 2 and = x\f + xla 
Jo 

where (j = 1)2) stands for the characteristic function of the set VFj in fl. It is 
easy to see that x'ji^) = Xji^l^) ^ ^p) where Xj is the characteristic function of 
the set Gj in . With the above notation, we may reformulate system (ini)-(ini) 
as a variational problem, that is, find it® G C^([0, T]; n C([0, T]; IJq( fl)^) 

and Pe G L^(0, T; L^(n|)) such that 

^{{u^t),p)) + (^xl^o'^u^t)+X 2 BoX/^{t)^ ■ Vpdx- J^xlPe 

^ lo (/j "^^ 2 -Si(t-'r)V^(r)^ ■Vipdx'^dT 

= m,ip)) ioT allipGH^in)^. 

The choice of a global test function (p G HQ{n)^ instead of two independent test 
functions on the two subdomains allows us to incorporate the interface conditions 
(frill-(IT31) into the weak formulation. Equation (12.11) requires G T^(0, T ; 
may be written under the equivalent divergence-free weak form 


dt^ 




^x!24|Vw^(t)+X^i?§V—(t) 

,du^ 


Vpdx 


x\A\{t - T)yu^T) F xlBlit - t)V—( r; 
/o Vn \ Di 

= ((r£,p)), all p G W, 


Vipdx j dr 
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or equivalently 

^ (p)) + {{Ae{t)u%t), (p))+ J - t)u%t), (fi)) dr = iiFe,(p)) V(/9 e W, 

( 2 . 2 ) 

Ae(t) and Be{t — r) (for 0 < t < t <T) being bounded operators defined from W 
into W by 

{{Ae{t)u^{t),^)) = +X 2 BoV^{t)j ■ Vifdx 

{{Be(t-T)u^{T),Lp)) = ^xfA®(t-T)VM®(r) +x^B®(t-T)V^(T)^ -Vtfdx. 

In order to ensure the existence and uniqueness of we use respectively El 
Theorem 3.4] and [HI Theorem 3.2] as follows. First, in view of [TTJ Theorem 3.4], 
the operator ^ +Ae is one-to-one from T^(0, T ; W) onto T^(0, T ; W'). Next since 
Beit — t) is bounded from W into IT', it readily follows from [HI Theorem 3.2] 
that (12.2p possesses a unique solution G L^(0,r;IT) with G L^(0,T;IT'). 
Therefore by setting and = X2^r> readily obtain a solution 

(Ue,Ve) G L^(0,T;\F) X L^(0,T;V)?) for the divergence-free formulation of (11.11) - 
(ini), and owing to the regularity of the data (see assumption (A3)), the following 
result is in order. 


Theorem 2.1. For each e > 0 and under Assumptions (Al)-(A3), there exists a 
unique triplet {ue,Ve^Pe) which possesses the regularity G IT^’°°(0, T; 14^), G 
L^iO,T;L^ni)^), Ve G tTi’“(0, T; iJ,) n 1.2(0, T; if), p, G L°°(0, T; 
and satisfies equations (fOl-dlTll. Moreover, G C\[0,T]-, L^{ni)^), G 
C([0,T];ile) and the following estimates hold: 

sup ^||M£(t)|f2(f2p -|- ||VM£(t)||^2(f2e) -|- ||r’£(t)||^2(f2E)^ ^ C, ( 2 - 3 ) 


ess sup 

0<t<T 



2 

L2(n=) 




sup 

0<t<T 


dUe 

dt 


it) 


< C and 


iVu, 


L^ni) 


lL2(n=) dt < C, 


l|PellL2(n|x(0,T)) — tA 

where C is a positive constant not depending upon e. 


(2.4) 

(2.5) 

( 2 . 6 ) 


Proof. It just remains to prove the existence and uniqueness of the pressure, to¬ 
gether with the uniform estimates. But the existence and uniqueness of p£ is derived 
exactly as in m Section 3.2] (see especially Theorem 3.4 therein). 

It remains to verify inequalities (I2.3I) - (I2.6I) . Multiply Eq. p.ip by (t) = u'^it), 
Eq. (11.21) by v^it) and sum up, then 


\Pi~d ^’ at ) 


(P' 


dvs 
2 dt 


{BfiS/Ue + BI* VUe, Vt!e)^ 2 (n: 


) 


-I A]; * Vltg, u'fij 
= (Plf T'^'e)L^{Ql) + iP29^'^e)L‘^(ni) 


(2.7) 
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where (, )l 2 stands for the inner product in But ()2.7p is equivalent to 


ld_ 

2 dt 

1 d 




L=(no 






2 dt * Vile, Vn£)^2,'Qej + {Bq'VVs Bf ^ V'Ue, 

= iPlf,Kit))LHQl) + iP29,Veit))^2(^ne) . 

Integrate (|2.8D with respect to t, 


( 2 . 8 ) 




L2(n=) 




L2(ni) 


+ * VMe, VMe)^ 2 (qs) + 2 /p (SgVWe + * Vtle, VtlE)^ 2 (Qs) dr 


= 2/o (pf/, <(i))L2(ne) + (p|£/,^^e(i))L2(n= 


dr. 


But 


/o 


(pf /, Me(t))i2(ne) + (p^g, We(i))L2(n=) 


< 




L2(nf 

‘ '"(pD’«'(^)"' 


+ 




z,2(n= 


dr 

dr 


L2(r2f) 


(p^)2t,^(T) 


L2(n|) 


dr. 


Making use of (Al) together with (I2.8I) - (I2.9|) we obtain 


(2.9) 


(pf) = <W 


L^ni) 


+ 


(p|)2Vg(t) 


L2(n|) 


■a||VMe(t)||i2(ne) +2a/Q ||Vi;e(T)||^2(n|) dr 


L2(n| 


djT (-^1 * 
2 

L2(Of) 


<Io( (Pi) = /L) ,2,0=^+ (P2)^9{t) 

\ i (f^l) . 

- 2 /o (^!* Vt;e,V^^e)i2(oe)dr+ /g + ||1^2(01)) 

where here the convolution * is taken with respect to t. But using Young’s inequal¬ 
ity, 

{A\ * VMe, VMe)^ 2 (Qe) 

= (fo ^i(-> ^ - r)Vu,(T)dT, Vu,(t)^ 

= fo (/of 'r)VMe(T) ■ Vue{t)dx^ dr 

< fo (c Pi(t - 'r)VMe(T)||^2(oe) + # l|VMe(t)||^2(oe)) dr 

< f ||VM^(t)||^2(oe) + C/o ||Af(t - T)VMe(r)lli2(ne) dr 

< f ||VMe(t)|li2(oe) + C/o ||Vwe(T)||i2(ne) dr, 


hence 


(A^ * Vme, Vtte)^2(f2e) < 2 ll^'*^e('^)|lL2(Qp dr. (2.11) 
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Also 


= 2 /o (fj (fni ■ yv^{T)dxj ds^ dr 

^ 2 /o (C'/J' l|Sf(^-s)Vi;e(s)|li2(ne)ds + /J' § ||VvjT)||^,(fje) ds) dr 
fo l|Vt’e(T)||i2(n|) dr + C J* (/J" ||Vi;e(s)|li2(oe) dsj dr, 


< a 


hence 


< a 


((BJ * Vi;e)(T),Vi;e(T))^2(s^e)dT 

£ ||VijaT)|li 2 (n|) dr + cj^ ||Vi;,(s) 


( 2 . 12 ) 


lL2(n|)' 


ds dr. 


{piV^KiO 

2 

+ 

L^{ni) 

0 

2 

+ - 
L2(n|) 2 

<c + Sl[ 


+ 

L2(Of) 

(p|)2t>£(T) 


Putting together (12.111) and (12.1211 . it emerges from (I2.10|) that 

^ ||VMe(t)||^2(oe) +a/g ||Vt>e(T)||^2(n|)dT 
2 

L2(0|) 

AC l|VME(T)||i2(f;e) + C/q" 11V (s) 11 ^2 (n|) ds]dT. 

Appealing to Gronwall’s inequality, 

2 

L2(nf) ^‘(i‘ 2 ) ' .Jo 

for all 0 < t < T and all e > 0. We recall that in all the above inequalities, C is a 
positive constant independent of both e and t. We therefore infer from (A2) that 

ll'*^£(^)llL2(n=) + ll'*^e(0llL2(f2S) + ll^'*J'£(^)llz,2(f2S)+ J II V'*J£(j')llL2(f2|) du < C (2.13) 


(pD = «£(i) \\'^'^eiT)\\l2^f^s)dT <C 


Now, from 


Ue{t) = / u'^{r)d7 

Jo 


we get 


Il“e(0lli2(ne) < Wu'.ir) 
thus, making use of ()2.13l) it holds that 


lL2(ne) 


dr 


ll'*^£(0llL2(f2e) < C, 

C > 0 being independent of both t and e. All the estimates (I2.3p - (I2.5I) follow 
thereby. Finally, for (12.6p . we proceed exactly as in the proof of m Section 3.2] 
(see also [Q] Lemma 2.3]). This completes the proof. □ 


In order to deal with the compactness result, let us recall the definition of the 
auxiliary functions we have just used in the preceding result: 

We{t) = [ v^{s)ds, = xfWe + xl^e and = x\p\ + X 2 P 2 . 

^0 


































ACOUSTICS OF ALMOST PERIODIC POROUS MEDIA 


9 


With this in mind, we easily see that the vector function {ue,We,ps) satisfies the 
weak formulation (fTTI) . (ina - (l2J6ll and (fTa-dlTll where 

^ 2 -^^ - div + J - s)V-^{x,s)dsj +Vpe = p^g 

(2.14) 

div = 0 in n 2 T (2-15) 

Me = Mig on rf 2 X (0,T) (2-16) 

where t = ^2 ^ (0j2^)- In Eq- (12.141) . represents the fluid displacement. In 
Section 5 we shall deal with the asymptotic behaviour of the above reformulation 
of our problem. It follows at once from Theorem 12.11 that the above system pos¬ 
sesses a unique solution (me, w^,pg) such that G W'^’°°{0,T;H^)r\H^{0,T-V^). 
Moreover it holds that 

Proposition 2.1. The sequence (m®)e>o is relatively compact in H^{0,T; 

Proof. First, in view of interface conditions (12.161) and (11.51) we have Vm® = 
yfVME -I- Next, the outer boundary conditions (11.61) imply m^ = 0 on 

i9fl. Thus m'^ G Hq{Q)^). Besides, the inequalities in Theorem 12.11 

ensure that the sequence (M'^)e>o is bounded in the space 

W = G L\0,T;H^{nf) : ^ G L^{0,T -and ^ G L^{0,T -, 

and the latter space is compactly embedded in (0, T ; L‘^{n)^). The result follows 
thereby. □ 


3. Almost periodic functions 

Let B(R.^) denote the Banach space of all bounded continuous (real-valued) 
functions on endowed with the sup norm topology. 

A mapping u : —>■ R is called an almost periodic function if u G B{R^) 

and further the set of all its translates {«(• + a)}ogRN has a compact closure in 
B{R^). We denote by AP(E.^) the set of all continuous almost periodic functions 
on . AP(]R'^) is a commutative C*-algebra with identity. An argument due to 
Bohr [5] (see also [3]) specifies that a function u G AP(]R^) if and only if u may 
be uniformly approximated by finite linear combinations of functions in the set 
{cos(fc • y), sin(fc ■ y) : k € R^}. It emerges from the above equivalent definition 
that AP(R^) consists of uniformly continuous functions. Furthermore, it is now 
well-known that AP(R^) enjoys the following properties: 

(P)i u{- -I- a) G AP(R^) whenever u G AP(M^) and for every a G 

(P )2 For each u G AP(R'^) the closed convex hull of {M(--l-a)}aGR'^ S(R'^) con¬ 

tains a unique complex constant M(u) called the mean value of u, and which 
satisfies the property that the sequence (u^)e>o (where u^{x) = u{x/e), 
X G R^) weakly ^-converges in L°°(R^) to M{u) as £ ^ 0. Moreover the 
following properties of M are in order: 

(i) M is nonnegative, i.e., M{u) > 0 for any u G AP(R^) with u > 0; 

(ii) M is continuous on AP(R^); 

(hi) M(l) = 1; 

(iv) M is translation invariant, i.e., M (u(--l-a)) = M{u) for all u G AP(R'^) 
and all a G R^. 
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We infer from the above properties that AP(]R^) is an algebra with mean value 
on [13]. We denote by /C its spectrum and by Q the Gelfand transformation 
on AP(R^). We recall that /C is the set of all nonzero multiplicative linear forms 
on AP(]R^), and G is the mapping of AP(]R'^) into C{IC) such that Giu){s) = s{u) 
(s G /C, u G AP(]R^)). The image G{u) of u will very often be denoted by u. 
We endow 1C with the relative weak* topology on (AP(]R^))' (the topological dual 
of AP(]R^)). In this topology, G is an isometric ^-isomorphism of the C*-algebra 
AP(K^) onto the C*-algebra C{JC). 1C being topologized as above, the mapping 
j : ^ /C given by j{y) = Sy (the Dirac mass at y) is continuous with dense 

range. 

In view of the properties (i)-(iii) in (P) 2 , the mean value admits an integral 
representation with respect to some Radon measure m (3 (of total mass 1) as 
follows: 

M{u) = [ G{u)d(3 for all u G AP(R^). 

Jk 

The following important result is worth recalling. 

Theorem 3.1. The topological space K, can he provided with a group operation 
under which it is a compact topological Abelian group additively written whose the 
Haar measure is precisely the Radon measure /3. Moreover, the mapping j : R^ —^ 
1C given by j{y) = Sy is a continuous group homomorphism satisfying the property 

rffu = Tj(^y-ju, all u G AP(R^) and all y G R^ (3.1) 

where^ denotes the Gelfand transformation on AP(R'^), TyU = u{-+y) andTj(^yyu = 
u{-+5y). 

Proof. The fact that 1C can be provided with a group operation is classically known; 
see e.g. |8]. For the proof of (EU, see the proof of Proposition 1 in [27]. □ 

In the sequel, 1C is viewed as a compact Abelian group additively written, 
equipped with the normalized Haar measure /3. Hence, as it is customary, we 
shall often write ds for dp. 

Remark 3.1. The group operation defined on 1C is actually the continuous exten¬ 
sion of the mapping 1C x 1C ^ 1C defined by (Sx, Sy) i—>■ Sx+y Since it is additively 
written, we are justified in writing Sx + Sy := Sx+y for x, y G R^. Hence, denoting 
by —s the symmetrization of s G we have in particular —Sy := S-y. 

Next, we introduce the space AP°°(M^) = {u G AP(R^) : DyU G AP(R^) for 
every a = (oi,..., ojv) G N^} where ■ 

Now, let H^p(R^) (1 < p < oo) denote the space of Besicovitch almost periodic 
functions on M^, that is the closure of AP(R^) with respect to the Besicovitch 
seminorm ^ 

||u|| = (limsup-— / \u{y)f dy] 

\r^+cx \.Dr\ JBr J 

where Br is the open ball of R'^ of radius r centered at the origin. It is known that 
is a complete seminormed vector space verifying i?^p(R^) C H^p(M^) 
for I < p < q < oo. Using this last property one may naturally define the space 
H^p(R^) as follows: 

B^p{R^) = {f eni<p+^BlpiR^)-. sup ||/||p<oo}. 

l<p<oo 
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We endow B^p{R^) with the seminorm [f]^ = supi<p<(^ ll/llpi which makes it a 
complete seminormed space. We recall that the spaces i3^p(R'^) (1 < p < oo) are 
not Frechet spaces since they are not separated. The following properties are worth 
noticing [191 Subsection 2.2] (see in particular Theorem 2.6 therein); see also [2^ 
Subsection 2.2]: 

(1) The Gelfand transformation Q : AP(]R^) —>■ C{IC) extends by continuity to 

a unique continuous linear mapping, still denoted by Q, of into 

LP(/C), which in turn induces an isometric isomorphism Qi, of B^p{M.^)/Af = 

onto LP(/C) (where M = {u € B^p(R^) : Q{u) = 0}). More¬ 
over if M e i?^p(M'^) n L°°(K'^) then Q{u) G L°°(/C) and < 

(2) The mean value M, defined on AP(IR'^), extends by continuity to a positive 
continuous linear form (still denoted by M) on B^p(IR.^) satisfying M{u) = 
fi^G(u)d/3 and M{u{- + a)) = M{u) for each u G i3^p(R^) and all a G 

Moreover for u G i?^p(M^) we have ||m||p = [M{\uf)]^^^. 

Spaces of almost periodic functions with values in a Banach space are defined 
in a natural way, we refer to [3] for details. Keep the following notations in mind: 
AP(R^;IR) = AP(R^) and B^p(M^;M) = B^p(K^). 

Let X Rt- denote the space R^ x R with generic variables (y, t). 

It holds that AP(R^+^) = AP(Rt-; AP(R^)) is the closure in 6(R^+^) of the 
tensor product AP(R^) 0 AP(Ri-) [121 Proposition 2.3]. To avoid heaviness of 
notation, we may sometimes set Ay = AP(R^), Ar = AP(Rt-) and A = AP(R^+^). 
Correspondingly, we will denote the mean value on A^ (^ = y, r) by M^. 

Now let 1 < p < oo and consider the 7V-parameter group of isometries {T{y) : 
y G M^} defined by 

T{y) : 6^p(R'^) ^ S^p(R'^), T(y){u+M) = u{-+ y) + M for u G Blp(M.^). 

Since AP(R'^) consists of uniformly continuous functions, {T{y) : y G R^} is a 
strongly continuous group in the following sense: T{y){u+N) u-\-M in i3^p(R^) 

as Ip] —0. Using the isometric isomorphism Qi we associated to {T{y) : y G R^} 
the following Wparameter group {T{y) : y G M'^} defined by 

T(p) : LP{IC) LP{IC) 

T{y)g,{u+M) = gi{T{y){u+M)) = 5i(w(- + y) + Af) for u G Blp{R^). 

The group {T{y) : y G R^} is also strongly continuous. The infinitesimal generator 
of T{y) (resp. r(y)) along the ith coordinate direction, denoted by d/dyi (resp. 
di) is defined by 

^ = lims^o {T{sei)u - u) in S^p(R^) 

(resp. diV = lims_j.o (T{sei)v — v) in LP{IC)), 

where we have used the same letter u to denote the equivalence class of an element 
u G B^p{]K^) in S^p(R^), = {Sij)i<j<N (Sij being the Kronecker (5). We 

collect here below some properties and spaces attached to d/dyi and di. We refer 
the reader to e.g. [H Subsection 3.1] for details. First, we denote by g the canonical 
mapping of i?^p(R^) onto S^p(R^), that is, g{u) = u + Af ioi u G U^p(R^). 

Let 1 < i < N. Then 
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(1) If u e AP(R^) is such that ^ e AP(]R^), then 


= e 



(3.2) 


(2) Gi o -^ = di o Qi provided that these mappings make sense. 

We define the Sobolev type space S^p(K'^) and its smooth counter-part I^ap(R^) 
as follows: 


= {u e : — e S^p(R^) VI < z < TV} 

^Ui 

and 

I?ap(R^) = p(AP°°(IR^)). 

First we have that X>ap(R^) is dense in ;B^p(]R^), 1 < p < oo. Next, endowed 
with the norm 


I l.p 


= hli: + 


N 

E 

i=l 


du 


dpi 


i/p 


(u e 


B^p(]R^) is a Banach space admitting 2 ?ap(R'^) as a dense subspace. We still write 
u either for G{u) if u G i?^p(IR^) or for Giiu) if u G S^p(IR^). 

We now define the appropriate space of correctors. Prior to that, we set 

= {u G : M{u) = 0}. 


We endow it with the norm 




du 




i/p 


(uGB 




Under this norm B^p(IR'^)/]R is unfortunately not complete. We denote by S^^p 
its completion with respect to the above norm and by J the canonical embedding 
of ,B^p(]R^)/]R into ,B^^p(M^). It can be easily checked that Pap(R'^)/R = {m G 
2?ap(R'^) : M{u) = 0} is dense in ,B^p(]R^)/]R. The following hold true: 


(Pi) The gradient operator Vy = 

extends by continuity to a unique mapping denoted by Vy = 
^#AP 


{ill; .47) ■ bIpC*")/® ^ (BSp(R")) 

(4l) 


1<2<N 


and 


N 


II#,P 


HE 




^ with 

the 

d 


d 

dpi 


dpi 



i/p 

du 

\ 


dpi 


for u G S^Ap( 




(P 2 ) The space J(X>ap(I^'^)/C) is dense in S^^p(]R^). 
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(P3) The mapping Vy is an isometric embedding of l3^^p{R^) onto a closed sub¬ 
space of so that is a reflexive Banach space. By 

duality we define the divergence operator divy : —>■ 

{p' =P/{P-^)) by 

(divyM,n) = — (u,yyv'j for v € 6^^p(K^) and u = (ui) S (S^p(K^))^, 
where (^u, VyV^ = Ylf=i Ijc Uidivdp. 

We may also define the Laplacian operator on B^p (denoted here by Ay) 

as follows: 

(^AyW,v) = (diVy(Vy?n),z;) = — (^VyW,S7yv'j for all v G B^^p(M^). 

Then it is immediate that Ayg{u) = Q{Ayu) for all u G AP°“(]R^) where Ay stands 
for the usual Laplacian operator on Ry . 

Let us finally give a word of warning. Since J is an embedding, this allows us 
to view S^p(]R^)/R (and hence 'Dap{R^)/R) as a dense subspace of B^^p{M.^). 
Thus we shall henceforth omit J in the notation if it is understood from the context 
and there is no risk of confusion. Accordingly, Vy and ^ shall still be denoted by 

Vy and ^ respectively. 


4. The S-convergence 

The letter E will throughout denote any ordinary sequence (c„)„gN with 0 < 
£„ < 1 such that e„ —>■ 0 as n —>■ 00 . We shall denote ^ 0 as n —>■ 00 ” merely 
by ”£ —>■ 0”. In what follows, we use the same notation as in the preceding section. 
We assume in this subsection that is any open subset of . 

Definition 4.1. Let 1 < p < 00 . (1) A sequence (u£)e>o C LP{Q) is said to weakly 
"E-converge in LP{n) to some Uq G B^p(R^)) if as £ —>■ 0, we have 

/ Us{x)f(x, — ]dx^ // uo{x, s)f{x, s)dxds 

Jn ^ JJcixic 

for every / G L^'{ri; AP{R^)) (1/p' = 1 — 1/p), where Uq = Qio Uq and f = Q o f. 
We express this by writing Uq in (PI)-weak S. 

(2) The sequence (ug)g^Q C U‘(Pl) is said to strongly E-converge in LP(PI) to 
some Mo € LP{Pl;B^p{R^)) if it is weakly E-convergent towards uq and further 
satisfies the following condition: 

ll^ellLp(n) ^ ll'*^o|lLP(nxK:) • 

We denote this by writing —>■ mq in L^(ll)-strong E. 

The following results are in order; see e.g. m Theorems 3.1 and 3.5] for the 
proof. 

Theorem 4.1. Let 1 < p < 00 . Let (ue)eeE be o- bounded sequence in L'p(PI). Then 
there exists a subsequence E' from E such that the sequence {uf)e^E' is weakly E- 
convergent in L'p(PI). 
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Theorem 4.2. Let 1 < p < oo. Let {ue)eGE be a bounded sequence in 
Then there exist a subsequence E' of E and a couple of functions 

{uo,ui) G X LPin;B]^lpiR^)) 

such that, as E' 5 e ^ 0, 


dus 

dxi 


Me —)• Mo in W"^'P{^l)-weak; 

w LP{n)-weak S, 1 < * < 

uXi dyi 


The next result deals with the S-convergence of a product of sequences. 


Theorem 4.3 ([^ Theorem 6]). Let 1 < p,q < oo and r > 1 be such that 
1/r = l/p+ 1/q < 1. Assume {ue)eeE C L^{fl) is weakly E-convergent in L‘^{fl) 
to some Mo G L‘^(fl;Bfp(M.^)), and {ve)e^E C LP{yi) is strongly E-convergent in 
LP{Ll) to some vo G LP{n;B^p{M.^)). Then the sequence {ueVe) e^E is weakly E- 
convergent in L'^{fl) to uqVq. 


As a consequence of the above theorem the following holds. 

Corollary 4.1. Let {ue)eGE C LP{VI) and (ve)eeE C L^ (12) fl L°^{Tl) (1 < p < oo 
and p' = p/(p — 1)) be two sequences such that: 

(i) Me — >■ Mo in LP{fl)-weak E; 

(ii) Me — )• Vo in L^ (LI)-strong E; 

(iii) (Me)eGE is hounded in L°^{Ll). 

Then Me Me —>■ mqMo in LP{Vt)-weak E. 


We are now able to deal with E-convergence of a convolution product of se¬ 
quences. Before we can do that, let us first and foremost define the convolution of 
functions defined on K,. Let p, g, m > 1 be real numbers satisfying ^ ^ = 1 + 

For M G LP{IC) and v G we define the convolution product u*v as follows: 

{u*v){s) = / u{r)v{s — r)dr, a.e. s £ 1C. 

JK 

Then * is well-defined and we have that u*v G L'"(/C) with the following Young 
inequality: 

\\^*' p \\ l ^{ k :) — II“IIlp(/c) ll^llL9(y) ■ 

Now let M G LP; LP(K.)) and v G L‘^(/C)). We define the two-scale 

convolution u * *v as follows: 


(m * =i=M)(a;, s) = / [{u{t,-)*v{x-t,-)) {s)]dt 


= u{t,r) * v{x — t, s — r)drdt, a.e. {x, s) £R^ X K. 

JR^ JK 

Then ** is well defined as an element of L'^fR.^ x K.) and satisfies 


11'“ * *^IIl'"(R"xK) - ll“llLP(R'^xy) ll^llL‘'(R"xy) ■ 

It is to be noted that if m G Lp{C l; Lp(IC)) where fl is an open subset of R^, and 
V G L‘^{R^; L'^{IC)), we may still define m * *m by viewing u as dehned in the whole 
of R^ X 1C; it suffices to take the extension by zero of u outside Cl x 1C. 
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Finally, for u G and v G S^p(R^)) we define the two- 

scale convolution denoted by ** as follows: u**v is that element of L"*(]R^; S™p(M^)) 
defined by 

Qi(u * *v) = u * *v. 

It also satisfies Young inequality. 

Now, let t G be fixed. Then {St/s) E>o is a net in the compact topological 
group /C, hence it possesses a weak=i= cluster point r G 1C. In the sequel we shall 
consider a subnet still denoted by {6t/E)e>o (if there is no danger of confusion) that 
converges to r in the topology of /C, i.e. 

(5 1 —>• r in /C-weak * as e ^ 0. (4.1) 

Let (ue)e>o be a sequence in L'p{CI) (I < p < oo) which is weakly E-convergent 
to uo G S^p(R^))- We will see that a macro-translation of Ug induces both 

micro- and macro-translations on the limit uq. This is the main goal of the following 
result whose proof can be found in [27] . 

Theorem 4.4 ([23 Proposition 3]). Let {ue)e>o be a sequence in LP{CV) (1 < p < 
oo) and let the sequence (ue)e>o be defined by 

Ve{x) = Ue{x + t) (x G n — t) . 

Finally, let uq G LP{Cl;B^p(R^)) and assume that —>■ uq in L^(fl)-weak E. If 

holds true then 

Ve —t Vo in. L^iCl — t)-weak E (4-2) 

where vq G LP{CI — t, S^p(]R^)) is defined by vq{x, s) = uq{x + t,s + r) for {x, s) G 
{Cl — t) X /C. 

We have now in hands all the ingredients to prove the convergence result involv¬ 
ing convolution product of sequences. 

Let p, g,TO > 1 be real numbers such that ^ ^ = 1 + ^- Let (u£)e>o C LP{Ct) 

and {ve)e>o C be two sequences. One may view as defined in the whole 

by taking its extension by zero outside Cl. Define 

{uE*Vg){x)= [ UE{t)vg{x - t)dt (a;GM^), 

Js.n 

which lies in L™{M.^) and satisfies the Young’s inequality 

||We * U£||^m(Q) < l|■^^e|lLP(0) lke|lL9(n) • (4-3) 

Here is that result. 

Theorem 4.5 (' [321 Theorem 2.6]). Let {ue)e>o and {ve)e>o be as above. Assume 
that, as e —>■ 0, Ue ^ uq in LP{Cl)-weak E and Ve —>• vq in L'^(E.^)-strong E, where 
Uq and Vq are in LP{CI; B^p{R^)) and L'^(M^; S^p(K^)) respectively. Then, as 
€ —^ 0, 

Ue * Ve ^ Uq * *vo in L'^{Cl)-weak E. 

Remark 4.1. In the case where the open set Cl is bounded and further p = 2 and 
q = 1, we have provided a full proof of this result in m (see especially the proof 
of Theorem 2 therein). 
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In this work, we shall use the evolving version of E-convergence involving space 
and time variable. Before stating it, however we need some preliminaries. 

Let Ay = AP(IR^) and Ar = AP(]Rt-). We know that A = AP(IR^+^) is the 
closure in i5(IR^+^) of the tensor product Ay 0 At- We denote by K.y (resp. /Ct, 
/C) the spectrum of Ay (resp. Ar, A). The same letter Q will denote the Gelfand 
transformation on Ay, Ar and A, as well. Points in ICy (resp. ICr) are denoted by 
s (resp. So)- The Haar measure on the compact group ICy (resp. K-r) is denoted by 
ds (resp. dso). We have 1C = ICy x ICr (Cartesian product) and the Haar measure 
on 1C is precisely the product measure (3 = ds® dso; the last equality follows in an 
obvious way by the density of Ay ® Ar in A and by the Fubini’s theorem. 

This being so, a sequence {ue)e>o C LP{Q) (1 < p < oo) is said to weakly 
E-converge in LP{Q) to some uq € L^((3; fi^p(R^+^)) if as e —>■ 0, we have 

/ uAx,t)f {x,t,—,-] dxdt ^ // uo{x,t,s,so)f{x,t,s,so)dxdtdsdso 

Jq \ JJqxk 

for every / G L^’'((5; AP(]R^+^)) (1/p' = 1 — 1/p). 

Remark 4.2. The conclusions of Theorems Q] and 14.21 are still valid mutatis 
mutandis in the present context (change there Cl into Q in Theorem 14.11 IF^’^(H) 
into LP{0,T-W^’P{n)) and W^'P{CI) x LP(H; S^Pp(K^)) into LP{Q,T-,W^^p{CI)) x 
LP(Q;S^p(M.;H^Pp«)))). 


Now, let T G M, and let {ue)^^Q C Lp (Q) (1 < p < oo) be a weakly E-convergent 
sequence in LP [Q) to uq G LP{Q-, S^p(]R^+^)). Set 

Ve{x, t) = Ue{x, t + t) for {x, t) G Q — {0, t) = Cl X (—r, T — t). 

Then Vg vq in Lp(Q — (0,r))-weak E where vq G Lp(Q — (0,r);S^p(K^+^)) is 
defined by 

vo(x, t, s, So) = uo{x, t + T,s,so + To), {x, t, s. So) G [Q - (0, r)] x 1C, 

the micro-translation xq being a weak=i: cluster point of the sequence {Sr/e)e>o- ^ 
similar conclusion holds for Theorem 14.51 mutatis mutandis. 

Consider a sequence (ue)e>o C LP{Q) (1 < p < oo) which is weakly E-convergent 
in LP{Q) to some vq G LP{Q-,B^p{R^~^^)). Define 

We{x,t) = / Ve{x,T)dT a.e. {x,t) G Q. 

Jo 

Then the following important consequence of Theorem 14.51 holds. 


Corollary 4.2. We have 

We —t Wo in LP{Q)-weak E as e —>■ 0 
where wq G Lp(Q; l3^p(R-^)) is defined by 

wo{x,t,s) = 

Proof. By definition, We{x,t) = (1 * Ve){x,t) where here, 1 stands for the constant 
sequence assuming value 1, and the convolution is taken with respect to the variable 


vo{x, T, s, ro)drodT, {x, t, s) G Q x ICy 


0 JKr 
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t. Theorem 14.51 tells us that Wg —>■ 1 * *vo in LP(Q)-weak S as e ^ 0 where by 
definition, 

(1 * *vo)(x,t,s,so) = / / vo{x,T,s,so -ro)drodT {{x,t,s,so) € Q X 1C) 

Jo JK-r 

which does not depend on sq as Vo(x, t, s, Sq — rQ)drQ = vq(x,t, s,rQ)drQ; 
remind that /Cr is a topological group. The proof is complete. □ 

Remark 4.3. It follows immediately that the limit of the sequence (we)e>o defined 
above does not depend on the microscopic temporal variable r. 

5. Asymptotic behaviour 

Let 0 be the characteristic function of the set S in S being considered in 
Subsection 1.1. We aim at studying the asymptotics of the sequence {u^,Vg,pg)g^Q 
defined in Section 2, under the additional assumption 
(A4) 

9 is almost periodic (see Subsection 1.1), (5.1) 

Ao,Bo G {Blp{R^))^\ Ai,Bi G AP{R^+Y^ and m G AP(R'^) (j = 1,2). 

(5.2) 

However, for some obvious reasons given in Section 2, we shall rather deal with 
the sequence {u^,w^,pg)^^Q satisfying (11.11) . (I2.14I) - (I2.16I) . (I1.5I) - (I1.7I) . 

Throughout this section we do not distinguish between a function u G i3^p(K"*) 
and its equivalence class u + A/” in S^p(]R’"). This will of course not lead to any 
confusion. Here H denotes a bounded open subset of with Lipschitz boundary 
dCl. We still set Q = Cl x (0,T), and the other notations are as in the preceding 
sections. 

5.1. Preliminary results. To begin with, it is important to recall that Xj (j = 
1, 2) is the characteristic function of the set Gj (defined in Section 1) in . With 
this in mind, the following result is the point of departure of the whole homoge¬ 
nization process. 

Lemma 5.1. Under assumption (15.11) we have 

Xj G Blp{R^) with M{xj) > 0 for j = 1, 2. (5.3) 

Proof. By virtue of m Corollary 3.2] (see also [3TJ Proposition 4.1 and Corollary 
4.2]) we have ^ ^ap(®^) -^(Xi) > 0- On the other hand, since X 1 +X 2 = 1 

a.e. in R.^, it follows at once that X 2 C B\p(R^) with M{x 2 ) > 0. But by the 
definition of the geometry of the domain, we have M{xi) < 1; nnd from the equality 
M{xi) + M{x 2 ) = 1, it readily follows that M(x 2 ) >0. □ 

The next result relies on the previous one and its easy proof can be found in |18] 
(see also M)- 

Lemma 5.2. Let {us)e>o be a sequence in L^{Q) (1 < p < 00 ) which weakly S- 
converges in LP{Q) towards uq G LP{Q-,B\p{R^~^^)). For j = 1,2, let Gj be the 
set defined in Section 1. Then, as e —> 0, 

■WeX^ uoXj L^(Q)-weak S. 
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Let E be an ordinary sequence defined as at the beginning of Section 4. As in 
Section 2, let 

=X\Ue + X2We, (5.4) 

Ug and ide be as in Section 2. Then is the displacement of the e-fluid-solid 
structure. We know from Section 2 that G iL^(0, T; iL,)with d'^u^/dt^ £ 
L^(0,T;L^(n)'^) = L'^{Q)^. Also, as in Section 2, we set = XiPi + X 2 P 2 (the 
density of the e-fluid-solid structure). By virtue of Proposition 12.11 we derive the 
existence of a subsequence E' of E and of a function Ug G E^(0,T; such 

that 

Ug in iL^(0,T;L^(n)^) strongly as i?' 9 £ —)• 0. (5.5) 

Invoking the diagonal process, we extract from {u^)g^E' another subsequence (still 
denoted by {u^)^^e') such that 

Ug in iL^(0, T; iLg (fl)^)-weak and -5> L^((5)^-weak. (5.6) 

It follows that Ug G with d'^ugjdt^ G Lp‘{Q)^. We infer from 

Theorem 14.21 and Remark 14.21 that (up to a subsequence of E' not relabeled) 

^ Vug + VyUi in L^{Q)^ -weak S (5.7) 

where Ui G L^(Q; SAp(®r; S^ap(®^))^)- Applying Lemma [521 we readily see 
that, as A' 9 £ ^ 0, 

xfViiE = xfVit® Xi(Vmo -I- Vyiti) in L^{Q)^ -weak E, (5.8) 


and 

= x^^w^ —>■ X2(^“o + Vyiti) in L^{Q)^ -weak E. 
We also observe that 

^ ^ ^ ™ L^(Q)^'-weak E 


(5.9) 

(5.10) 


when E' 9 £ ^ 0, where Mr{vi) = dui/dt in G 2 . Indeed for the last convergence 
result above, we have, in view of (15.6p . 


X2' 


- dw 


dt ^X2^inL^(Q)^^-weakE. 


Also, there exists vi G such that 

^ ™ L^((5)^"-weak E. 

We seek the relationship between Ui and Vi in G 2 . First, invoking Corollary 14.21 

XzVtUe = X 2 y ^ MriVyVi)^ {T])dr) in -weak E. 

The uniqueness of the limit implies 


X 2 (Vmo + VyMi) = X 2 J + MriVyVi)^ (r])dr]. 

Taking the weak derivative with respect to t in the above equality. 


i9iti 

'W 


X2Vy^ = X2^r(VyUi), Or X2Vy 


i9iti 

'W 


= 0 . 
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The algebra of continuous almost periodic functions being ergodic, it follows from 
[34l Lemma 2.2] that — Mr{vi) agrees in G 2 with a constant in K^, hence 

- Mrivi) = 0 in G 2 (5-11) 

since the sole function ip in 6^p(Rr; satisfying Vyip = 0 is zero. It 

comes that iti G iL^(0, T; 6 ap(®t; S^ap(®^))’^)) depend on r in 

G 2 . This will be of great interest in the next subsection. Indeed, this will allow us 
to select test functions not depending on r in G 2 . 


5.2. Passage to the limit. In this section we pass to the limit in the variational 
formulation of (11.11) . (I2.14I) - (I2.16L (ll.5l) - (ll.7D (which is equivalent to the one of 
dnn-dni)). For the sake of simplicity, we shall omit throughout this subsection 
to precise that G' 9 e —5> 0 when dealing with a convergence result. Bearing 
this in mind, let = (ipQ,'ipi) be such that ip^ = {ipQ)i<k<N S C^iQ)^ and 
'^1 = ii’i,k)i<k<N G (Cg°(Q) (g) with dxpi/dr = 0 in G 2 (that is, 

ipi does not depend on r in G 2 ). We define = i/’o + by 

^e(x, t) = ■0o(a;, t) + £%pi t, j, for (x, t) G Q. 

Then we have fFg G C“((5)^ and, taking \I/g as a test function in the variational 
formulation of dtn), (Eni-diiiD, (iT3])-(fT:7i). we get (after using the interface 
conditions (II. 5p and (I2.16|) l 

- Jq • ^edxdt + Jq Xi + Af* Vme) • V^^dxdt 

+ SQXl{Biy^+Bl*y^)-V^,dxdt (5.12) 


- !q XlPe div ^edxdt = fgiXlPlf + xlplde) ■ ^edxdt. 

Our main objective is to pass to the limit in (15.121) . Let us consider the terms in 
(I5.12|) separately. We first deal with the term fg p^^g^ ■ '^edxdt. We have 


d‘^u^ 

• ^edxdt = - 


>Q 


'Q 


dt dt 


-dxdt 


-du’^ 


dtPo 

dt 


d-iPi 

dr 


+ £ 


dipi 

dt 


dxdt 


and 


Hence 


d'it, dlpo dipi . ^2tn,\N 


dt 




dt 


+ ^ in L^(Q)^-strong S. 


dr 


. 92 ? 


dt^ 


"^edxdt 


p- 


^duQ j dipQ dip I 


QxK 


dt 


dt 


dr 


dxdtd/3 


where p = XiPi + X 2 P 2 - But ■ ^^dpi = 0 (as Uq does not depend on r, the 

microscopic temporal variable) and Hqxk'p^W ' ^^dxdtdp) = pds\ /g ■ 
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d'tpf 

dt 


■dxdt = — (^J/c 7>ds^ Jq • 'ipgdxdt. Hence, on letting Pq = 'pds, we have, 




'Q 


d'^uo 

dt^ 


■ 'tpQdxdt. 


Now, dealing with the second term of the left-hand side of (I5.12L one easily shows 
that 

Vi/>Q -f '^y'4>i in L^{Q)^ -strong E. 


V’®', 


(5.13) 


Putting (15.131) and (15.81) together, we get by Corollary 14.II that 

xfVite • V’®'e -5> Xi (Vmq -I- VyMi) • (Vi/Jq + in L^((5)-weak E. 

Hence, using Aq as test function (recall that Aq G ((H^p(R^) (~l so 

that by [25l Proposition 8], it is an admissible test function in the sense of [25l 
Definition 5]) leads to 


XiAqVus ■ VAigdxdt 


AqUu ■ Ih'i’dxdtdp 


QxK 


where, setting u = (ug, Ui), we define Dm = (D_;m)i<j<jv with 1 


T>ju'^)l<k<N 


andDjM'^ = -^+djUi {dju\ = Gi [dui/dyj)) with Ug = {ug)i<k<N (and the same 
definition for Mi). Still for the same term, let us now deal with the part involving 
convolution. First we know that A^ —>■ Ai in LP{Q)^ -strong E (1 < p < oo). 
Hence by virtue of Theorem 14.51 we obtain 

Al * Xi^'*^e Ai * *Xi (Vmq -I- VyMi) in -weak E 

(just take p = 1 above). Therefore by repeating the reasoning above, we quickly 
arrive at 


QxK 


Xii^i * *Dm) • DAfdxdtdp. 


/ Xl(^l * ^“e) ■ 

JQ 

Thus 

/ Xl + ^1 * ^“e) - 

J Q J JQxK 

Arguing as for the preceding term and accounting of (15.111) and of the fact that 
does not depend on r in G 2 , we get 

,dwp_ „dw 

1h 


'Q 




dt 


•V®', 


Xi (AqDm-I-A l = 1 = *Dm) • D®'da;(i<(i/3. 


dxdt ^ ff x2iBo^^u+Bi**^I]iu)-]D>"^dxdtdf3. 
JJQxK at 


As regards the terms with the pressure, one has 

/ X2Pe div = / x^Pediv-i/jQfiidt-I- / X2Pe(divj/■0i)®da;(it(5.14) 

JQ JQ JQ 

-\-e / X 2 Pe(div' 0 i)^da;dt. 

JQ 

Owing to (12.61) . we consider p € L‘^{Q; Hip(R^+i)) such that 

xlPe X 2 P in T^(Q)-weak E 
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up to a subsequence of E' not relabeled. Set po{x, t) = X 2 P{^^ ■Sj so)d/3 for a.e. 

{x, t) G Q. Then passing to the limit in (j5.14|) yields 


JQ 

It also holds that 


X^Pe div "if^dxdt po div i/igdxdt + 


X2P‘i^^4>idxdtdl3. 


QxK 


/ {x\p\f + X2P29) ■ ^edxdt ^ (XiPif + X2P29) ■ ipodxdtdp. 

JQ JJQxK 

Finally, from the equation div = 0 in x (0, T) which is equivalent to 


Xl div ^ = 0 in Q, we get 


QxK.^ 


^ ( duQ -p- dui \ ^ , 

Yn div —-h div—^ 'ibdxdtds = 0. 

\ ot at 


(5.15) 


Indeed, since ui does not depend on r in G 2 , we choose a test function ip G 
Cpp^iQ) (8> AP“(IR^), and we have 


/» Q'uij 

J div -^dxdt = 0. 


'Q 

Passing to the limit above and using ()5.11l) . we are led to (I5.15|) . 

We have just verified that the triplet {uo,Ui,p) determined earlier solves the 
variational problem 

Po Iq ■ ^odxdt + Xi{j^o^u + Ai * *Du) • ID'i/dxdtd/d 

+ //qxk X20o§i^u + Bi* • 'D'ifdxdtd/3 

- /g Po div ipgdxdt - ff^,^ X2Pd^^^idxdtdf3 = J^F ■ xp^dxdt, 

//qxk:„ X 2 (div ^ + dA^) ^dxdtds = 0 

for all (i/^o.V’i) e C^iyQ f X [Co"”(g) 0 AP°“(K^+i)]^ and i/' G 


Co°“(Q)®AP= 


where 


Po = 


/ iXiPi + X2P2)ds and F = {XiPif+ X 2 P 29 )ds. 

Jk„ Jk„ 


Taking successively {xpQ.xpi) = and (■0o,’0i) = (0, lAi) in (15.161) . we get 

the following equivalent system: 




dt^ 




—div XiiM'OyU + Ai = 1 = *Dyu) + X 2 Bo^ByU + Bi * *—ByU +Vpo = F, 






dt 




dt 


\\ 


(5.17) 


-divj; Xii^o^yU + Ai = 1 = =i=DyM) + X 2 BQ—OyU + * *—DyW +Vy(x2P) = 0, 


dr^ 

I a- '^“0 , dui 

^,|d,v —+d,v,— 


dt 


= 0 , 


dUn . 

Mo(a:,0) = -^{x,0) = 0 


(5.18) 

(5.19) 

(5.20) 


where here, DyM = (Dm) = Vuq + V^Mi. 

Lemma 5.3. The system (I5.17I) - (I5.20I) has a unique solution. 
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Proof. The proof follows the same lines of reasoning as the proof of [121 Lemma 
5], □ 


Our aim now is to derive the effective equations for Uq. In order to do so, we 
need to consider firstly (j5.18l) and (15.191) which will help us to express iti and p in 
terms of Uq. We shall seek {ui,p) under the form Ui = Un + Ui 2 where Un and 
(mi 2 ,p) solve respectively Eq. (15.2111 and (I5.22|) below: 

— divy (^Xii^o^yU + ^1 * =i=Dyi6)) = 0 (5.21) 


f -divj^ {x2iBo-§i^yU + Bi* + Vy{x2P) =0 .. 

1 x,(div^+dl^,^)=0. ^'-22) 

Eq. (I5.2ip allows us to find Ui in Gi while (15.221) permits us to determine Ui in 
G 2 and p in G 2 , so that itn = Xi^^i and M 12 = X 2 '*^i- Let us first deal with (15.211) . 
For ^ = {iij)i<i,j<N S be fixed, consider the cell equation 

( Find u(^) e such that 

I -dI^j,(klo(C + V,,u(a)+Ai**(C + Vj,u(C)))=OinGi (5.23) 

[ u{(f) = 0 in G 2 . 


In view of the assumptions on Aq and Ai, we can verify, using mi Theorem 6.4], 
that (I5.23P possesses a unique solution. Taking in (I5.23P the special ^ = Vito(x,t) 
{{x,t) G Q) and comparing the resulting equation with (I5.21|l . we get by the 
uniqueness of u{^) that ui{x,t) = u{VuQ{x,t)) in Gi, which define a mapping 
u{Vuq) : {x,t) u{\7uo{x,t)) from Q into 

Now, as for (15.221) . fix again ^ = {^ij)i<i,j<N in and consider this time the 
following equation 

f -divy (go(g + Vyu(O) + Bi* + ^yv{^))) + ^yTT{0 = 0 in G 2 24 i 
\ tr(^) + divyi;(^) = 0 in G 2 and v{^) = 0 in Gi 


where tr(^) = is the trace of Let B]^p(G 2 ) denote the strong closure 

in of the set V 2 = {'0 G (I?a(R^))^ : ^yW = 0 and i/; = 0 in Gi}. 

Then following the standard methods for solving compressible Stokes equations 
mi Chap. 1], it also follows from mi Theorem 6.4] that (I5.24p possesses a unique 
solution (u (^), TT (^)) £ B\p (K^ ; B^p (G 2 )) x B\p (K^ ; B\p (R^ ) /R) . This means that 
v{^) is unique up to an additive function f £ Sap(Rt ;Bap (G 2 )) such that VyC = 0 
in G 2 and 7r(^) is unique up to an additive function which is constant in G 2 . Now, 
fixing ^ = ^'S/uq{xA) {ix,t) G Q) in (15.241) and repeating the same arguments as 
above, we end up with 


Ui = V 


dt 


S/Ur 


in G 2 and p = tt 


dt 


Vut 


in G 2 


where the functions v (^Vmq) and tt (^Vito) are defined in the same manner as 
u (Vmq). This being so the functions Ui and p are therefore defined by 


wi = Xiw(Vtio) + X2V 


9^ 


and p = X27r 



We are now able to derive the effective equations for the equivalent medium. 
This is the goal assigned to the next subsection. 
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5.3. Homogenization result. We begin by defining the effective coefficients. Let 
the matrices Cq = {cQ)i<i,j<N and Ci = be defined as follows: for 


any f 


II 

[ Xi(^o(? + du{0) + Ai* + du{0))df3 

JK 

CiC = 

f MBo{^ + di^)) + Bi**{^ + d^)))d/3, 

JK 

and the function 


HO = [ X2^{0dfi- 

JK 


Then by virtue of the uniqueness of u(^), n(^) and 7r(^) (for a given ^), Co, Ci and 
h are well defined. Now, substituting Ui by + X 2 '^ and p by 

X 2 '^ (^Vmo) in (15.171) and taking into account (15.201) . we get the model equation 
for our equivalent medium, viz. 

r - div (CoVmo + CiV^) + Vh (V^) =FmQ 

< tio = 0 on X (0,T) (5.25) 

i 0) “ “o(a:, 0) = 0 in n. 

The fact that (j5.25|) possesses a unique solution is an easy consequence of Lemma 

[El 

We are therefore led to the homogenization result. 

Theorem 5.1. Assume that (A1)-(A4) hold. For cache > 0, let {u^^v^^pi;) be the 
unique solution to (ICT-lfTTl). Let denote the global displacement field defined 
by (15.41) . Then the sequence (M'^)e>o strongly converges in to Uq 

where Uq is the unique solution to (|5.25l) . 

Proof. We just need to prove the convergence of the whole sequence (tt^)£>o. But 
this is an obvious consequence of the uniqueness of the solution to (15.251) . □ 
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